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Abstract 


Land subsidence is a significant hazard affecting many plains of Iran, including the 
Eyvanakey Plain in Semnan Province. The Eyvanakey Plain has a high potential for 
subsidence due to the influence of climate and topography. Given the importance of this issue, 
the present research aims to evaluate land subsidence in the Eyvanakey Plain and analyze the 
factors contributing to its occurrence. 

In this study, Sentinel-1 radar images, Landsat satellite imagery, and data from piezometric 
wells in the region were utilized as primary research data. The key tools employed include 
GMT, ArcGIS, ENVI, and SPSS software. The study was conducted in two main phases: 
first, an assessment of land subsidence in the region from 2016 to 2022, and second, an 
analysis of the impact of groundwater depletion and land use changes on subsidence. 

The results indicate that the study area experienced subsidence ranging from 28 to 533 mm 
over the six-year period, which is a significant amount. Spatial analysis of the subsidence 
reveals that the highest rates occurred in the southern parts of the Eyvanakey Plain. These 
areas also experienced an annual drop in groundwater levels of approximately 2 meters, 
highlighting groundwater depletion as one of the primary causes of subsidence. 
Furthermore, based on the results, the extent of man-made areas and agricultural lands 
increased by 6.9 km? and 2.7 km?, respectively, between 1992 and 2022. Given that the 
highest levels of subsidence were observed in these land use areas, the development of 
agricultural lands and man-made areas has been identified as the main cause of subsidence. 


Key words: Land Subsidence, Groundwater Depletion, Eyvanakey Plain, Land Use changes, 
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حکیده 


فرونشست زمین مخاطره‌ای است که بسیاری از دشت‌های ایران ازجمله دشت ایوانکی در استان سمنان را دربر گرفته 
است. دشت ایوانکی تحت تأثیر وضعیت اقلیمی و توپوگرافی» پتانسیل بالایی از نظر وقوع مخاطره فرونشست دارد. با 
توجه به اهمیت موضوع. در این پژوهش به برآورد وضعیت فرونشست دشت ایوانکی و تحلیل عوامل مؤثر در وقوع 
آن پرداخته شده است. در این تحقیق از تصاویر راداری سنتینل ۰۱ تصاویر ماهواره لندست و اطلاعات مربوط به چاه- 
های پیزومتری منطقه به‌عنوان مهم‌ترین داده‌های تحقیق استفاده شده است. مهم‌ترین ابزارهای پژوهش. GMT‏ 
ENVI ArcGIS‏ و SPSS‏ بوده است. با توجه به موضوع و اهداف مورد نظر» در این پژوهش ابتدا به ارزیابی 
وضعیت فرونشست منطقه در طی سال‌های ۲۰۱۶ تا ۲۰۲۲ و سپس به ارزیابی تأثیر افت منابع آب زیرزمینی و تغییرات 
کاربری اراضی در فرونشست رخ داده پرداخته شده است. بر اساس نتایج حاصله. محدوده مطالعاتی در طی دوره 
زمانی ۶ ساله بین YA‏ تا ۵۳۳ میلی‌متر فرونشست داشته که رقم قابل‌توجهی می‌باشد. آنالیز مکانی فرونشست رخ داده 
Sle‏ این است که بیش‌ترین میزان فرونشست مربوط به مناطق جنوبی دشت ایوانکی بوده است که این مناطق نیز 
سالانه: با افت سطح آب زیرزمینی حدود Y‏ متر مواجه شده است. بر این اساس» یکی از دلایل فرونشست رخ داده؛ 
افت منابع آب زیرزمینی بوده است. همچنین بر اساس نتایج حاصله. کاربری نواحی انسان‌ساخت و اراضی کشاورزی 


در طی سال‌های ۱۹۹۲ تا ۲ به ترتیب با ۶/۹ و ۷/۲ کیلومترمربع افزایش مواجه شده‌اند و با توجه به اینکه بیش - 
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ترین میزان فرونشست (بیش از ۵۰۰ میلی‌متر) در محدوده این کاربری‌ها بوده است. بنابراین توسعه اراضی کشاورزی 
و نواحی انسان‌سانعت. عامل اصلی وقوع فرونشست منطقه بوده است. 

کلیدواژه‌ها: فرونشست زمین.کاهش منابع آب زیرزمینی. دشت ایوانکی.تغیبرات کاربری اراضی.تصاویر راداری 
tes‏ تل کات 


۱- مقدمه 


امروزه روند افزایشی جمعیت و دخالت‌های نابجای OL!‏ در محیط زمینه را برای گسترش مخاطرات طبیعی 
فراهم آورده است )2019 (Ganjaeian,‏ ازجمله مخاطرات طبیعی که در طی سال‌های اخیر بسیاری از جوامع 
را درگیر کرده است» مخاطره فرونشست است )2009 (Delinom et al.,‏ بنا به تعریف سازمان زمین‌شناسی 
ایالات‌متحده. فرونشست عبارت است از فروریزش و يا نشست سطح زمین که به علت‌های متفاوتی در مقیاس 
So»‏ روی می‌دهد(2010 (Amighpey, Arabi & Talebi,‏ فرونشست زمین مسئله‌ای جهانی و پدیده‌ای 
مورفولوژیکی است(2011 (Hsieh et aL,‏ این پدیده متأثر از فعالیت‌های انسانی و عوامل طبیعی است که 
ممکن است مبدل به مخاطره و تهدیدی برای انسان و دستاوردهای انسانی شود( Saffari, Jafari & Tavakoli‏ 
(Sabour, 2018‏ فرونشست زمین وابسته به وضعیت زمین‌شناسی. اقلیمی و هیدروژئومورفولوژی مناطق است 
و عوامل انسانی مانند برداشت از منابع GT‏ زیرزمینی و سازه‌های انسانی نقش تشدیدکننده دارند( Li et al,‏ 


.(2023; Hasibuan et al., 2023 


مناطق مختلفی از جهان تحت تأثیر وضعیت طبیعی و عوامل انسانی» در معرض مخاطره فرونشست قرار 
دارند )2023 Malik, Kumar, Perissin & Pradhan, 2022; Jiang et al.,‏ وضعیت هیدرواقلیمی ایران سبب 
شده است تا بخش‌های زیادی از آن در معرض این مخاطره باشند. در واقع» با توجه به اینکه بخش زیادی از 
وسعت ایران را مناطق خشک و نیمه‌حشک دربرگرفته است و این مناطق با کمبود منابع آب سطحی مواجه 
هستند. بنابراین این مناطق در معرض فرونشست قرار دارند. ازجمله مناطقی که در معرض مخاطره 
فرونشست هستند» دشت‌های استان سمنان ازجمله دشت ایوانکی است. بررسی‌های اولیه صورت گرفته و 


مین گزارگی‌های آراته شده از سر cela llo‏ ملف GUS‏ داه است: که دشت ایوانکی با تشاطره 


۱۳۰۳ نشریه جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره ۲ پائیز‎ YA 


فرونشست مواجه است و اثرات ol‏ به‌صورت ترک‌خوردگی اراضی کشاورزی. ترک‌خوردگی سازه‌های 
انسانی و غیره نمایان شده است. با توجه به اهمیت موضوع و اثرات زیانبار فرونشست ej‏ است تا به 
بررسی مخاطره فرونشست در این دشت پرداخته شود که در این پژوهش به این مهم پرداخته شده است. 

در ارتباط با موضوع موردمطالعه تحقیقات مختلفی در سطح ایران و جهان صورت گرفته است که ازجمله 
آن‌ها می‌توان به کیم و همکاران )2007 (Kim, Kim, Kim, Won & Moon,‏ اشاره کرد که با استفاده از روش 
سری زمانی Olga PSI‏ فرونشست منطقه شهری بوسان را ۳۰ میلی‌متر در سال برآورد کرده‌اند. بزانو و 
همکاران )2015 (Bozzano et al.,‏ با استفاده از روش تداخل‌سنجی راداری» ميزان فرونشست دشت اکه E‏ 
ارا حدود ۸۰ میلی‌متر در سال برآورد کرده‌اند. چن و همکاران(2015 (Chen et al.,‏ با استفاده از روش سری 
زمانی PSI‏ به ارزیابی میزان فرونشست شهر پکن پرداختند. دا لیو و توسی )2018 (Da Lio & Tosi,‏ په 
ارزیابی فرونشست زمین در دشت ساحلی فریولی ایتالیا پرداختند. نتایج تحقیق آن‌ها نشان داده است که 
میزان فرونشست این دشت بیش از ۵ میلی‌متر در سال بوده است. بخارائی و همکاران( Bokhari et al.,‏ 
L (2023‏ استفاده از روش سری زمانی PSI‏ میزان فرونشست زمین در شهر گوادار" پاکستان را در طی دوره 
زمانی ۲ ساله بین ۳۲ تا ٩۳‏ میلی‌متر فرونشست برآورد کر ده‌اند. خان و همکاران( & Khan, Faiz, Gadea‏ 
Ahmad, 2023)‏ با cole!‏ از روش تداخل‌سنجی راداری» میزان فرونشست در حوضه پیشاور پاکستان را 
حدود ۳/۲۳ سانتی‌متر در سال برآورد کرده‌اند. در ایران نیز صفاری و همکاران(2018 (Saffari et al.,‏ با 
استفاده از روش سری زمانی SBAS‏ فرونشست NUS‏ دشت شهریار-کرج را در طی سال‌های ۲۰۰۳ تا 
۰ ۱۳۶ میلی‌متر برآورد کرده‌اند. بابایی و همکاران(2017 (Babaei, Khazaei & Qasere Mobarakeh,‏ با 
استفاده از روش‌های PS‏ و SBAS‏ نشان دادند که در قسمت جنوب غربی شهر تهران و در مناطق شهرداری 
۱٩ (A ۰۱۷ A‏ و ۲۱ فرونشست معناداری در فاصله زمانی این یک سال اتفاق افتاده است. اصغری 
سراسکانرود و محمدزاده شیشه گران(2021 (Asghari Saraskanroud & Mohamadzadeh Shishegaran,‏ با 


استفاده از روش تداخل سنجی راداری» ماکزیمم ميزان فرونشست زمین در oboe eu‏ حدود ۱۱/۵ سانتی - 


1 Acque Albule 


? Gwadar 


عباس‌زاده و همکاران» برآورد نرخ فرونشست دشت ایوانکی و تحلیل نقش فعالیت‌های انسانی در وقوع آن ۳۹ 


متر در سال برآورد کرده‌اند. قره‌جلو و همکاران Gharechelou, Akbari Ghoochani, Golian & Ganji,‏ 
)2021( با استفاده از روش تداخل‌سنجی راداری. ماکزیمم میزان فرونشست زمین در دشت مشهد را حدود 
۰ سانتی‌متر در سال برآورد کرده‌اند. اسدی و همکاران( Asadi, Ganjaeian, Javedani & Ghaderi Hasab,‏ 
1 به ارزیابی میزان فرونشست دشت ایوانکی با استفاده از تصاویر راداری پرداختند. نتایج این پژوهش 
نشان داده است که دشت (Sly!‏ در طی سال‌های ۶ تا ۲۰۱۸ ماکزیمم ۲۳ سانتی‌متر فرونشست داشته 
است. حسین‌زاده و همکاران )2023 (Hosseinzadeh, Akbari, Javanshiri & Mohammadpour,‏ با استفاده از 
روش تداخل‌سنجی راداری» میزان فرونشست زمین در دشت مرکزی شهرستان قائن را بین ۲/۶ تا VIA‏ 
سانتی‌متر در سال برآورد کرده‌اند. محرابی و همکاران(2023 (Mehrabi, Karimi & Khalesi,‏ با استفاده از 
روش تداخل‌سنجی راداری OUS‏ دادند که میزان فرونشست دشت جیرفت از ۱۱ سانتی‌متر در سال ۲۰۱۴ به 
۳ سانتی‌متر در سال ۲۰۲۲ افزایش ah‏ است. در راستای تحقیقات پیشین صورت گرفته. هدف از این 
پژوهش. برآورد وضعیت فرونشست دشت ایوانکی و تحلیل عوامل مؤثر در وقوع Of‏ است. در این تحقیق بر 
خلاف تحقیقات پیشین ازجمله اسدی و همکاران(2021 (Asadi et al.,‏ دوره زمانی ۶ ساله در نظر گرفته شده 
است و علاوه بر محاسبه نرخ فرونشست. ارتباط آن با افت منابع آب زیرزمینی و تغییرات کاربری اراضی نیز 
ارزیابی و تحلیل شده است. 

-Y‏ مواد و روش‌ها 
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ف o allia‏ ماش MS‏ وف اھ کے در شب اسان تماق cst E abd‏ سوه 
۰ کیلومترمربع وسعت. از نظر تقسیمات سیاسی در محدوده شهرستان گرمسار قرار دارد و مهم‌ترین شهر 
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۰ متری از سطح boo‏ قرار دارد. همچنین این دشت از نظر اقلیمی sim‏ مناطق خشک تا نیمه‌عشک 
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Fig.1. Location map of the study area 
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-Y-Y‏ روش انجام پژوهش 


در این تحقیق از تصاویر راداری سنتینل c‏ تصاویر ماهواره لندست و اطلاعات مربوط به جاه‌های پیزومتری 


منطقه Ul pea‏ مهم‌ترین داده‌های تحقیق استفاده شده است. مهم‌ترین ابزارهای پژوهش. GMT"‏ (به‌منظور تهیه 


و SPSS‏ (به‌منظور انجام محاسبات مورد نظر) بوده است. با توجه به موضوع و اهداف مورد o‏ این پژوهش 


در چند ale yo‏ انجام شده است که در ادامه به تشریح آن‌ها پرداخته شده است: 


+ Generic Mapping Tools 


۳۱ 
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مرحله اول (ارزیابی ميزان فرونشست منطقه): در این مرحله با استفاده از تصاویر سنتینل | مربوط به دوره 


زماتی ۶ساله ار ۲۰۱۶ نا ژانویه و BAS adi o ghey‏ فیزان tude Catia‏ 


ایوانکی محاسبه شده odd‏ در این مرحله از ۶۱ تصاویر راداری استفاده oJ‏ است (جدول تصاویر مورد 


استفاده دارای حالت مداری صعودی» نوع IW‏ فرمت SLC‏ و پولاریزاسیون ۷ هستند پس از تهیه تصاوین 


به‌منظور پردازش تصاویر از سیستم عامل لینوکس " و نرم‌افزار GMT‏ استفاده شده است که برای این منظور در 


تصاویر انجام شده است. همچنین در این مرحله آستانه همبستگی تصاویر ۰/۲ تعریف شده است. پس از 


انجام پیش‌پردازش‌های لازم بر che‏ بیس‌لاین " زمانی تصاویر نقشه‌های اینترفروگرام منطقه تهیه شده است. 


در این پژوهش ابتدا ميزان جابجایی سالانه محاسبه شده اشیت و درنهایت جابجایی A‏ صورت گرفته از 


سال ۲۰۱۶ تا ۲۰۲۲ محاسبه شده است. 


حدول \— تاریخ تصاو بر راداری مورد استفاده 


Table 1- Date of the radar images used 


ردیف تاریخ ردیف تاریخ ردیف تاریخ ردیف تاریخ 

Date Num. Date Num. Date Num. Date Num. 
2020.09.17 47 2019.03.27 32 2017.10.15 17 2016.01.06 1 
2020.10.23 48 2019.05.02 33 2017.11.20 18 2016.02.23 2 
2020.11.28 49 2019.06.07 34 2017.12.26 19 2016.03.18 3 
2021.01.03 50 2019.07.13 35 2018.01.19 20 2016.04.11 4 
2021.02.08 51 2019.08.18 36 2018.02.24 21 2016.05.29 5 
2021.03.16 22 2019.09.23 37 2018.04.01 22 2016.07.04 6 
2021.04.21 53 2019.10.29 38 2018.05.07 23 2016.09.02 7 
2021.05.27 54 2019.12.04 39 2018.06.12 24 2016.10.20 8 
2021.07.02 55 2020.01.09 40 2018.07.18 25 2016.12.07 9 
2021.08.07 56 2020.02.14 41 2018.08.23 26 2017.01.24 10 
2021.09.12 57 2020.03.21 42 2018.09.28 27 2017.03.01 11 

! Linux 


2 Baseline 
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2021.10.18 58 2020.04.26 43 2018.11.03 28 2017.04.06 12 


2021.11.23 59 2020.06.01 44 2018.12.09 29 2017.05.12 13 


2021.12.17 60 2020.07.07 45 2019.01.14 30 2017.06.29 14 


2020.01.10 61 2020.08.12 46 2019.02.19 31 2017.08.04 15 


2017.09.09 16 


مرحله دوم (تحلیل ارتباط فرونشست منطقه با افت منابع آب زیرزمینی): در این مرحله با استفاده از اطلاعات 
مربوط به وضعیت افت سطح آب چاه‌های دشت ایوانکی میزان تغییرات این چاه‌ها محاسبه شده و سپس 
نقشه وضعیت افت سطح آب زیرزمینی در منطقه موردمطالعه تهیه شده است. در این مرحله ابتدا اطلاعات 
مربوط به ۱۳ چاه پیزومتری در محدوده مطالعاتی تهیه شده و سپس با استفاده از نرم‌افزار SPSS‏ میزان افت 
سالانه هر چاه محاسبه شده است. پس از محاسبه میزان افت سالانه هر چاه با استفاده از روش درونیابی 
Spilne‏ (بر اساس مطالعات کتابخانه‌ای صورت گرفته. این روش تناسب بیش‌تری برای این موضوع دارد) و 
از طریق نرم‌افزار ArcGIS‏ نقشه افت سالانه آب زیرزمینی منطقه تهیه شده است. پس از تهیه نقشه میزان 
افت سالانه col‏ زیرزمینی منطقه با استفاده از ضریب همبستگی پیرسون, ارتباط افت منابع آب زیرزمینی با 


فرونشست منطقه ارزیابی شده است. 


مره سوم led)‏ ارتباط فرولخسنت als‏ تغیبرات eligi‏ در این مرحله با استفاده از pala‏ 
ماهر ازه clade (Pus) ad‏ کازبری Gls calif‏ در be‏ سال‌های ۶۱۹۹۲ ۲۲۲ هه شاه انیت 
o AT cd‏ وات کین Pediat snis ws dcos eo) EL Det SU‏ 
ایی می vatis al‏ کازبری اراق Yes ey AY. ual do us oai alte‏ و ix ag Yu‏ 
است. پس از تهیه تصاویر بیش‌پردازش‌های لازم شامل تصحیحات هندسی و تصحیحات رادیومتریک بر 
روی تضاویر انجام شده است و سپس با استفاده از روش حداکثر احتمال d‏ نقشه‌های کازبری اراضی منطقه 
تهیه شده است. پس از تهیه نقشه‌های کاربری اراضی و آنالیز تغییرات صورت گرفته, به ارزیابی وضعیت 
فرونشست زمین در هر کاربری و ارتباط تغییرات کاربری اراضی با ميزان فرونشست رخ داده پرداخته شده 


است. 


1 Maximum Likelihood 


yY 


عباس‌زاده و همکاران. برآورد نرخ فرونشست کس ایوانکی و تحلیل نقش فعالیت‌های انسانی در وقوع ol‏ 


۳۴ نشریه جغرافیا و مخاطرات محیطیء جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


حدول -Y‏ مشخصات تصاویر مورد استفاده 


Table 2- Specifications of the images used 


ردیف تاریخ ماهواره سنجنده 

Sensor Satellite Date Num. 
TM Landsat 5 1992 1 
TM Landsat 5 2002 2 
ETM Landsat 7 2012 3 
OLI Landsat 8 2022 4 


۳نتایج و بحث 
۳-۱"ارزیایی میزان فرونشست منطقه در طی سال‌های ۲۰۱۶ تا ۲۰۲۲ 
در این پژوهش به‌منظور ارزیابی فرونشست منطقه از تاریخ ۲۰۱۶/۰۱/۰۶ تا ۲۱۲۲/۰۱/۱۰ از ۶۱ تصویر 
راداری استفاده شده است و بر اساس بیس OR‏ زمانی» 2 ae‏ تصویر به‌منظور تهیه نقشه‌های اینترفرو گرام 
انتخاب شده است AY J te)‏ پس از انتخاب زوج تصاویر اینترفروگرام‌های مورد نظر در این بازه زمانی 
تهیه شده Cs‏ (شکل AY‏ 

جدول ۳- زوج تصاویر انتخاب شده جهت تهیه نقشه‌های اینترفروگرام 


Table 3- Selected image pairs for generating interferogram maps 


o las » شما‎ I 
: Image pair زوج تصاویر‎ = Image pair „gla زوج‎ 


Image pair زوج تصاویر‎ 
ER Num. Num. Num. 


2020.02.14 2020.03.21 41 2018.02.24 2018.04.01 21 2016.01.06 2016.02.23 1 
2020.03.21 2020.04.26 42 2018.04.01 2018.05.07 22 2016.02.23 2016.03.18 2 
2020.04.26 2020.06.01 43 2018.05.07 2018.06.12 23 2016.03.18 2016.04.11 3 
2020.06.01 2020.07.07 44 2018.06.12 2018.07.18 24 2016.04.11 2016.05.29 4 
2020.07.07 2020.08.12 45 2018.07.18 2018.08.23 25 2016.05.29 2016.07.04 5 
2020.08.12 2020.09.17 46 2018.08.23 2018.09.28 26 2016.07.04 2016.09.02 6 
2020.09.17 2020.10.23 47 2018.09.28 2018.11.00 27 2016.09.02 2016.10.20 7 
2020.10.23 2020.11.28 48 2018.11.03 2018.12.00 28 2016.10.20 2017.12.07 8 
2020.11.28 2021.01.03 49 2018.12.00 2019.01.14 29 2016.12.07 2017.01.24 9 
2020.11.28 2021.01.03 50 2019.01.14 2019.02.19 30 2017.01.24 2017.03.01 10 
2021.01.03 2021.02.08 51 2019.02.19 2019.03.27 31 2017.03.01 2017.04.06 11 
2021.02.08 2021.03.16 52 2019.03.27 2019.05.02 32 2017.04.06 2017.05.12 12 
2021.03.16 2021.04.21 53 2019.05.02 2019.06.07 33 2017.05.12 2017.06.29 13 
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2021.04.21 2021.05.27 54 2019.06.07 2019.07.13 | 34 2017.06.29 4 14 


2021.05.27 2021.07.02 55 2019.07.13 2019.08.18 | 35 2017.08.04 2017.09.09 15 


2021.08.07 2021.09.12 56 2019.08.18 2019.09.23 36 2017.09.09 2017.10.15 16 


2021.09.12 2021.10.18 57 2019.09.23 2019.10.23 37 2017.10.15 2017.11.20 17 


2021.10.18 2021.11.23 58 2019.10.29 2019.12.04 38 2017.11.20 2017.12.26 18 


2021.11.23 2021.12.17 59 2019.12.04 2020.01.00 39 2017.12.26 2018.01.19 19 


2021.12.17 2022.01.10 60 2020.01.09 2020.02.14 40 2018.01.19 2018.02.24 20 
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)۲۰۱۶/۰۵/۲۹ - ۲۰۱۶/۰۷/۰۴ نمونه‌ای از نقشه اینترفروگرام‌های منطقه (زوج تصویر‎ -۲ JSS 
Fig.2. An example of interferogram maps of the area (image pair 2016.05.29 — 2016.07.04) 
نقشه میزان فرونشست‎ SBAS پس از تهیه نقشه‌های اینترفرو گرام منطقه. با استفاده از روش سری زمانی‎ 
شده اسٿ. بر اساس نتایج حاصله محد‌وده مطالعاتی در‎ ag Ye YY/* MY G (۱۲ ۱ ole ME منطفقه از تاریخ‎ 


Ub‏ دوره زمانی ۶ ساله بین YA‏ تا ۵۳۳ میلی‌متر فرونشست داشته که رقم قابل‌توجهی است. بر این اساس 


می‌توان گفت که دشت ایوانکی به‌طور متوسط سالانه حدود ۸ سانتی‌متر فرونشست دارد که این مسئله بیانگر 


بحرانی بودن و در ver‏ مخاطره بوده این ua qus‏ همچنین نتایج بررسی 4 C‏ پراکنش مکانی 


۳۶ نشریه جغرافیا و مخاطرات محیطیء جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


فرونشست رخ داده Sly‏ این است که بیش‌ترین میزان فرونشست مربوط به مناطق جنوبی دشت ایوانکی و 


بیش‌تر در حدفاصل روستاهای جنت‌آباد و چشمه نادری بوده است (شکل AP‏ 
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--. جاده اصلی‎ Main Road 28 فرونشست (میلیمتر)‎ 533 ; 9 


شکل ۳- نقشه تجمعی فرونشست منطقه از تاریخ ۲۰۱۶/۰۱/۰۶ تا ۲۰۲۲/۰۱/۱۰ 


Fig.3. Cumulative subsidence map of the area from the date 2016.01.06 to 2022.01.10 


۳-۲-بررسی وضعیت افت منابع آب زبرزمینی منطقه موردمطالعه 


یکی از مهم‌ترین عوامل وقوع فرونشست. افت منابع آب زیرزمینی است. در واقع» در ue‏ سال‌های اخیر 
روند افزایشی جمعیت و استفاده بیش‌ازحد از منابع آب زیرزمینی سبب افت سطح منابع آب زیرزمینی شده و 
همین مسئله زمینه را برای وقوع فرونشست در بسیاری از دشت‌های کشور ازجمله دشت ایوانکی فراهم 
آورده است. در این پژوهش به‌منظور بررسی وضعیت CH‏ سطح آب‌های زیرزمینی دشت ایوانکی از 
اطلاعات سازمان آب و منطقه‌ای استان سمنان استفاده شده است که برای این منظور میزان افت سطح آب ۱۳ 
چاه محاسبه شده است (جدول ۴). بر اساس نتایج حاصله. بیش‌ترین میانگین افت سالانه آب زیرزمینی با 


۱/۹۵ متر مربوط به جاه رج حیدر (جاه شماره ^( بوده است. بعدازاین جاه نیز جاه‌های تما مزرعه 


عباس‌زاده و همکاران» برآورد نرخ فرونشست دشت ایوانکی و تحلیل نقش فعالیت‌های انسانی در وقوع آن ۴۷ 


بهشتآباد ole)‏ شماره ۱۳) و چاه باغ دانش ole)‏ شماره ۷) به ترتیب با میانگین افت سالانه ۱/۷۵ و ۱/۵۱ 
متر» دارای بیش‌ترین میزان افت آب زیرزمینی بوده‌اند. همچنین کم‌ترین میانگین افت سالانه آب زیرزمینی با 
۱ متر مربوط به چاه چنداب ole)‏ شماره ۱) بوده است. در این پژوهش پس از محاسبه میانگین افت 
YL.‏ چاه‌های منطقه موردمطالعه. به‌منظور تهیه نقشه میزان افت سطح آب در منطقه موردمطالعه اطلاعات 
مربوطه وارد نرم‌افزار ArcGIS‏ شده و سپس با استفاده از روش درونیابی Spline‏ نقشه میزان افت سالائه سطح 
آب زیرزمینی در این دشت تهیه شده است (شکل ۴) که بر اساس Ol‏ پیش ترین میزان افت آب زیرزمینی 
مربوط به مناطق جنوبی دشت ایوانکی بوده است. 
جدول ۴- موقعیت و میانگین افت سالانه چاه‌های منطقه موردمطالعه 


Table-4 Location and average annual decline of wells in the study area 


موقعیت در سطح آب در سال سطح آب در سال میانگین افت سالانه 
نقشه موقعیت محلی دوره زمانی مبنا ۱۳۹۸ (متر) 
Average Water level in Water levelin Time period ^ Local location Location‏ 
annual decline the year 2019 ۰ the base year on the map‏ 
ندات 
P 1‏ 1993-2019 120 144 0.91 
Chandab‏ 
باغ انا 
us 2‏ 1993-2019 61 89 1.06 
Bagh Anari‏ 
Bened 3‏ 1993-2019 72 100 1.08 
Fuel pump‏ 
پارک بخشداری 
i 4‏ 1993-2019 142 170 1.07 
Park‏ 
آبا 
e 5‏ 1993-2019 118 136 1.14 
Jahad Abad‏ 
سیالک 
6 : 1993-2019 43 71 1.09 
Sialk‏ 
باغ Sala‏ 
Er 7‏ 1993-2019 96 122 151 
Garden‏ 
Ms 8‏ 1993-2019 91 142 1.95 
Borj Heydar‏ 
جنتآبا 
bin 9‏ 1993-2019 32 64 123 
Jannat Abad‏ 
شرق چشمه نادی 
East of 10‏ 1993-2019 6 42 14 


Cheshme Nadi 


۱۴۰۳ شماره ۳ پائیز‎ Y نشریه جغرافیا و مخاطرات محیطی» جلد‎ YA 


جنوب چشمه نادی 
South of Nadi 11‏ 1993-2019 4 38 131 
Cheshme‏ 
جنوب بهشت اباد 
South of 12‏ 1993-2019 4 37 1.26 
Behesht Abad‏ 
شمال مزرعه 


1.57 78 51 2002-2019 بت اباد‎ 14 
Northth of 
Behesht Abad 
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LL محدوده شهری‎ Urban Boundary 


شکل ۴- نقشه درونیایی شده افت آب زیرزمینی در منطقه موردمطالعه 
Fig.4. The interpolated map of groundwater decline in the study area‏ 
پس از نقشه درونیابی افت منابع آب زیرزمینی» به‌منظور ارزیابی ارتباط بین میزان فرونشست منطقه با افت 
سطح آب چاه‌های موردمطالعه. همبستگی بین آن‌ها محاسبه شده است که نتایج یافته‌ها Sly‏ ارتباط قوی و 


مثبت UA)‏ بین bol‏ می‌باشد (شکل ۵). با توجه به نتایج به دست آمده. یکی از عوامل اصلی در وقوع 


فرونشست منطقه افت منابع آب زیرزمینی بوده تیاه 


عباس‌زاده و همکاران برآورد نرخ فرونشست دشت ایوانکی و تحلیل نقش فعالیت‌های انسانی در وقوع آن ۳۹ 


Correlation همبستگی‎ 


y = 308.29 +2 
R? = 0.679 


فرونشست 
9 
5 


s) 


r3 
8 


Subsidence(mm) (> 


0 0.5 1 1.5 2 


افت آب زیرزمینی (متر) Groundwater Decline(m)‏ 
شکل ۵- نمودار وضعیت همبستگی بین افت سطح آب چاه‌ها با میزان فرونشست زمین 
Fig.5. Graph of correlation between groundwater level decline and subsidence rate‏ 


۳-۳-ارزیابی روند تغییرات کاربری اراضی در منطقه موردمطالعه 


با توجه به اینکه یکی از عوامل de‏ در وقوع فرونشست. تغییرات کاربری اراضی است. در این بخش به 
بررسی روند تغییرات کاربری اراضی محدود مطالعاتی و ارتباط آن با فرونشست منطقه پرداخته شده است. 
در این پژوهش به‌منظور بررسی تغییرات کاربری اراضی محدوده مطالعاتی از تصاویر ماهواره لندست مربوط 
به سال‌های ۸۱۹۹۲ ۲۰۰۲ ۲۰۱۲ و ۲۰۲۲ استفاده شده است. پس از تهیه تصاویر مورد نظر و انجام 
تصحیحات و پیش‌پردازش‌های لازم. نقشه کاربری اراضی محدوده مطالعاتی در طی سال‌های ۱۹۹۲ تا ۲۰۲۲ 
تهیه شده است (شکل ۶). بر اساس نقشه‌های تهیه شده, بخش زیادی از وسعت محدوده را مراتع و سپس 


۵۰ نشریه جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 
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شکل ۶- نقشه کاربری اراضی منطقه در طی سال‌های ۱۹۹۲ تا ۲۰۲۲ 


Fig.6. Land use map of the area from 1992 to 2‏ 
بررسی نقشه‌های کاربری اراضی محدوده مطالعاتی نشان داده است که این محدوده در طی سال‌های اخیر با 
تغییرات زیادی مواجه شده است. بر اساس نتایج حاصله. کاربری نواحی انسان‌ساخت در سال 1۹۹۲ ۱/۷ 
کیلومترمربع وسعت داشته است که این میزان در طی سال‌های ۰۲۰۰۲ ۲۰۱۲ و ۲۰۲۲ به ترتیب به YIN‏ ۷/۵ و 
۶ کیلومترمربم افزایش ail‏ است که این مسئله بیانگر روند توسعه این کاربری در طی سال‌های pel‏ بوده 
است. کاربری اراضی کشاورزی در سال ۱۹٩۲‏ ۴۱/۳ کیلومترمربع وسعت داشته است که این میزان در طی 
سال‌های ۲۰۰۲ ۲۰۱۲ و ۲۰۲۲ به ترتیب به ۴۴/۷ ۴۵/۲ و ۴۸/۵ کیلومترمربع افزایش یافته است که این 
مسئله بیانگر این است که کاربری اراضی کشاورزی نیز در طی سال‌های pH!‏ دارای روند افزایشی بوده است. 
کاربری مراتع و اراضی pl‏ در سال ۱۹۹۲ ۳۸۶/۲ کیلومترمربع وسعت داشته است که این میزان در طی 
سال‌های ۲۰۰۲ ۲۰۱۲ و ۲۰۲۲ به ترتیب به ۳۸۰/۷ ۳۷۶/۵ و ۳۷۲/۱ کیلومترمربع کاهش يافته است (جدول 
۵ و شکل ۲). با توجه به نتایج حاصله» کاربری نواحی انسان‌ساخت و اراضی کشاورزی در طی سال‌های 
۲ تا ۲۰۲۲ به ترتیب با ۶/۹ و ۷/۲ کیلومترمربع افزایش یافته مواجه شده‌اند و کاربری مراتع و اراضی بایر 


L‏ \/¥\ کیلومترمربع کاهش مواجه شده است. مجموع نتایج حاصله از این بخش نشان داده است که تحت 
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gla با‎ Eel eel کاس و‎ CONI ES ITE PE 
اشنا همچنین روند توسعه اراضی‎ à EA افزایش فشار بر زمین و تشدید فرونشست در منطقه‎ ¿eL مسئله‎ 
کشاورزی نیز سبب افزایش میزان بهره‌برداری از منابع آب زیرزمینی و درنتیجه وقوع فرونشست در دشت‎ 

TS‏ شده است.. در شکل ۸ تصویر از اثرات فرونشست در دشت Sly!‏ نشان داده شده است. 
حدول ۵- مساحت کاربری‌های اراضی محد و ده مطالعاتی در ub‏ سال‌های ۲۳ تا ۲۰۲۲ برحسب کیلومتر مربع 


Table 5-The area of Land use types within the study area from 1992 to 2022 (km?) 


۶ کا 
cu»)‏ س دربری 1992 2002 2012 2022 
Land use Type Row‏ 
Ks‏ انسان‌ساخحت 
petes a gaa 1‏ 17 3.8 7.5 8.6 
Human-Made‏ 
2 کشاورزی Agricultural‏ 413 44.7 45.2 48.5 
مراتع و اراضی p‏ 
Sk 05 3862 Rangelands and barren 3‏ 
lands‏ 
مراتع a‏ کشاورزی a‏ انسان iod Bess‏ 
Human-Made Agricultural Rangelands‏ 400 
350 
300 3 
250 1 
۷ 
3j 200‏ 
150 
z 100‏ 
P B‏ 
o .HB EN —B =.‏ 
2022 2012 2002 1992 
Year Juu‏ 


شکل ۷- نمودار مساحت کاربری‌های اراضی محدوده مطالعاتی در طی سال‌های ۱۹۹۲ تا ۲۰۲۲ 


Fig.7. Graph of the area of land use types within the study area from 1992 to 2022 
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شکل ۸- تصویر از اثرات فرونشست در دشت ایوانکی 


Fig.8. A picture of the effects of subsidence in the Eyvanakey Plain 
مجموع نتایج حاصله از این بخش نشان داده است که تحت تأثیر روند افزایشی جمعیت» روند توسعه‎ 
که این مسئله باعث افزایش فشار بر زمین و تشدید‎ abl فیزیکی نواحی سکونتگاهی و انسان‌ساخت افزایش‎ 
فرونشست در منطقه شده است. همچنین روند توسعه اراضی کشاورزی نیز سبب افزایش میزان بهره‌برداری‎ 
از منابع آب زیرزمینی و درنتیجه وقوع فرونشست در دشت ایوانکی شده است. در این پژوهش به‌منظور‎ 
ارزیابی ارتباط میزان فرونشست رخ داده در محدوده مطالعاتی با نوع کاربری‌های اراضی. وضعیت فرونشست‎ 
اطلاعاتی هر کاربری (کاربری‎ AY رخ داده در هر کاربری محاسبه شده است. به‌منظور انجام این کار ابتدا‎ 
کاربری (بر اساس نقشه تجمعی فرونشست منطقه در طی‎ Of سال ۲۰۲۲) جدا شده و سپس نقش فرونشست‎ 
سال‌های ۲۰۱۶ تا ۲۰۲۲) استخراج شده است. بررسی وضعیت میزان فرونشست رخ داده در کاربری مراتع و‎ 
۴۷۹ تا‎ YA داده است که میزان فرونشست این کاربری در طی سال‌های ۲۰۱۶ تا ۲۰۲۲ بین‎ GLE ub اراضی‎ 
میلی‌متر بوده است. بررسی وضعیت میزان فرونشست رخ داده در کاربری نواحی انسان‌ساخت نشان داده‎ 
است که میزان فرونشست این کاربری در طی سال‌های ۲۰۱۶ تا ۲۰۲۲ بین ۵۱ تا ۵۰۳ میلیمتر بوده است.‎ 
همچنین نتایج ارزیابی میزان فرونشست رخ داده در کاربری اراضی کشاورزی نشان داده است که میزان‎ 
پر آین‎ ۵ ASD) eal بردہ‎ din cla 0۵۳۳ EPA ین‎ Y Gc Y YA کازیری فرط سال‌های‎ ul فروتکست‎ 
اساس می‌توان گفت مناطقی که دارای کاربری اراضی کشاورزی هستند» دارای بیش‌ترین میزان فرونشست‎ 


هستند و بنابراین بین نوع کاربری کشاورزی E‏ فرونشست دح cools‏ ارتباط مستقیمی وجود دارد؛ بنابراین 


oY 
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طور مستقیم سبب افزایش بهره‌برداری از منابع آب زیرزمینی و در نتیجه وقوع مخاطره فرونشست شده است. 
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شکا 4— نقشه ميزان 3j‏ دسست کاربری اراضی کشاورزی در ub‏ سال‌های ۶ تا ۳۰۳۲ 


Fig.9. Map of land subsidence rates in agricultural land use from 2016 to 2 


۴-نتیجه گیری 


موقعیت قرا رگیری دشت اپوانکن سبب شده است تا این دشت پتانسیل فرونشست بالایی داشته باشد. در این 


سال‌های ۶ YYY LG‏ ارزیابی یله heal‏ بر اساس نتایج حاصله محد‌وده مطالعاتی TI.‏ دوره زمانی ۶ 


دشت ایوانکی به‌طور متوسط سالانه حدود ۸ سانتی‌متر فرونشست دارد که این مسئله Sly‏ بحرانی بودن و 


در معرض مخاطره بوده این دشت است. همچنین نتایج بررسی رضعیت پراکنش مکانی فرونشست ej‏ داده 


۵۴ نشریه جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


بیانگر این است که بیش ترین میزان فرونشست مربوط به مناطق جنوبی دشت ایوانکی و بیش‌تر در حدفاصل 
روستاهای جنت‌آباد و چشمه نادری بوده است. در این پژوهش به‌منظور بررسی عوامل موثر در فرونشست 
رخ داده» به بررسی وضعیت افت منابع آب زیرزمینی و تغییرات کاربری اراضی در دشت ایوانکی پرداخته 
شده است. بر اساس نتایج حاصله. میانگین سالانه افت آب چاه‌های پیزومتری منطقه بین ۰/۹۱ تا ۱/۹۵ متر 
بوده است و بر اساس نقشه تهیه شده» بیش‌ترین میزان افت سطح آب زیرزمینی, مربوط به مناطق جنوبی 
دشت ایوانکی بوده است. نتایج همبستگی بین میزان فرونشست منطقه با افت سطح آب چاه‌های موردمطالعه 
بیانگر ارتباط بین قوی و مثبت (۰/۶۸) بین آن‌ها می‌باشد. بر این اساس» همانند دشت کرج-شهریار ( Saffari‏ 
(et aL, 2018‏ و دشت مشهد(2021 «(Gharechelou et al.,‏ یکی از عوامل el‏ در وقوع فرونشست منطقه. 
افت منابع آب زیرزمینی بوده است. همچنین نتایج تغییرات کاربری اراضی OUS‏ داده است که کاربری نواحی 
انسان‌ساخت و اراضی کشاورزی در Gb‏ سال‌های ۱۹۹۲ تا ۲۰۲۲ به ترتیب با ۶/۹ و ۷/۲ کیلومترمربع افزایش 
يافته مواجه شده‌اند و کاربری مراتع و اراضی بایر با ۱۴/۱ کیلومترمربع کاهش مواجه شده است. با توجه به 
اینکه بیش‌ترین میزان فرونشست مربوط به کاربری اراضی کشاورزی و نواحی انسان‌ساحت بوده است؛ 
بنابراین می‌توان گفت که همانند دشت همدان(2023 «(Ganjaeian, Asadi, Menbari & Ebrahimi,‏ تحت 
تأثیر روند افزایشی جمعیت» روند توسعه فیزیکی نواحی سکونتگاهی و انسان‌ساخت افزایش «3b‏ که این 
مسثله باعث افزایش فشار بر زمین و تشدید فرونشست در منطقه شده است. همچنین روند توسعه اراضی 
کشاورزی نیز سبب افزایش میزان بهره‌برداری از منابع آب زیرزمینی و در نتیجه وقوع فرونشست در دشت 


ایوانکی شده است. 
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